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An exact study of charge-spin separation, pairing fluctuations and pseudogaps is carried out by 
combining the analytical eigenvalues of the four-site Hubbard clusters with the grand canonical and 
canonical ensemble approaches in a multidimensional parameter space of temperature (T), magnetic 
field (h), on-site interaction (U) and chemical potential (n). Our results, near the average number of 
electrons (A^) ~ 3, strongly suggest the existence of a critical parameter Uc{T) for the localization 
of electrons and a particle-hole binding (positive) gap A''~'^{T) > at f/ > Uc(T), with a zero 
temperature quantum critical point, Uc{0) = 4.584. For U < Uc{T), particle-particle pair binding is 
found with a (positive) pairing gap A^(r) > 0. The ground state degeneracy is lifted at f/ > Uc(T) 
and the cluster becomes a Mott-Hubbard like insulator due to the presence of energy gaps at all 
(allowed) integer numbers (1 < A < 8) of electrons. In contrast, for U < Uc(T), we find an electron 
pair binding instability at finite temperature near (A) « 3, which manifests a possible pairing 
mechanism, a precursor to superconductivity in small clusters. Rigorous criteria for the existence of 
many-body Mott-Hubbard like particle-hole and particle-particle pairings, spin-spin pairing, (spin) 
pseudogap and (spin) antiferromagnetic critical crossover temperatures, at which the corresponding 
pseudogaps disappear, are also formulated. In particular, the resulting phase diagram consisting of 
charge and spin pseudogaps, antiferromagnetic correlations, hole pairing with competing hole-rich 
((A) — 2), hole-poor ((A) = 4) and magnetic ((A) = 3) regions in the ensemble of clusters near 1/8 
filling closely resembles the phase diagrams and inhomogeneous phase separation recently found in 
the family of doped high Tc cuprates. 

PACS numbers: 65.80.-(-n, 73.22.-f, 71.10.Fd, 71.27.-(-a, 71.30.+h, 74.20.Mn 



I. INTRODUCTION 



Understanding the effects of electron correlation and 
pseudogap phenomena in doped oxides, 

including the cuprate superconductors is one of the most 
challenging problems in condensed matter physics 7]. 
Although the experimental determination of various in- 
homogeneous phases in cuprates is still somewhat con- 
troversial jgj, the underdoped high Tc superconductors 
(HTSCs) are often characterized by crossover tempera- 
tures below which excitation pseudogaps in the normal- 
state are seen to develop 0. In these materials, the 
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spectral weight begins to be strongly suppressed below 
some characteristic temperature Tg that is higher than 
the superconducting crossover temperature Tc. There 
are many experiments supporting a highly nonuniform 
hole distribution leading to the formation of hole-rich 
and hole-poor re gion s in doped La2-a:Sra;Cu04 and other 
cuprate HTSCs . This electronic phase separation is 
expected to be mostly pronounced at low hole concentra- 
tions. Recently, strong experimental evidence has been 
found for 'electronic phase separation' in La-cuprates 
near optimal doping into separate, magnetic and super- 
conducting phases 

The relevance of the Hubbard model for studies of the 
HTSCs has been the focus of intensive research and de- 
bated for quite some time with no firm conclusions up 
to now. Even though the small Hubbard clusters do 
not have the full capacity to describe the complexity of 



copper ions and their ancillary oxygens detected in the 
HTSC materials, it is still argued that this model can 
capture the essential physics of the HTSCs T|. However, 
beyond one and infinite dimensions, there is no exact so- 
lution currently available for the Hubbard Hamiltonian. 
It is also known that in the optimally doped cuprates, the 
correlation length of dynamical spin fluctuations is very 
small ,12], which points to the local character of elec- 
tron interactions in the cuprates. Therefore, exact mi- 
croscopic studies of pairing, crossover and pseudogaps, by 
using analytical diagonalization of small Hubbard clus- 
ters, which account accurately for short-range dynamical 
correlations, are relevant and useful with regard to under- 
standing the physics of the HTSCs. Our exact analytical 
solution appears to be providing useful insight into the 
physical origin of the high energy insulator-metal and low 
energy antiferromagnetic crossovers, electron pairing and 
spin density fluctuations in the superconducting phase. 

The following questions are central to our study: (i) 
When treated exactly, what essential features can the 
simple Hubbard clusters capture, that are in common 
with the HTSCs? (ii) Using simple cluster studies, 
is it possible to obtain a mesoscopic understanding of 
electron- hole/elcctron-electron pairing and identify vari- 
ous possible phases and crossover temperatures? iii) Do 
these small clusters (coupled to a particle bath) contain 
important features that are similar to large clusters and 
thermodynamical systems? 

Our work has uncovered important answers to the 
questions raised above. This is a follow-up to our re- 
cent study, in which rigorous criteria were found for the 
existence of microscopic quantum critical points (QCP), 
Mott-Hubbard (MH) and Neel type bifurcations, and cor- 
responding critical temperatures of crossovers and var- 
ious phases in finite-size systems |l3j . Small 2- and 
4-site clusters with short range electronic correlations 
provide (unique) insight into the thermodynamics and 
exact many body physics, difficult to obtain from ap- 
proximate methods. In particular, we show that these 
small Hubbard clusters, in the absence of long range or- 
der, exhibit particle-particle, Mott-Hubbard like particle- 
hole or antiferromagnetic-paramagnetic instabilities in 
the ground state and at finite temperatures. 

In addition, the 4-site (square) cluster is the basic 
building block of the Cu02 planes in the HTSCs and 
it can be used as a block reference to build up larger 
superblocks in 2D of desirable sizes by applying Clus- 
ter Perturbation Theory (CPT) Ti], non-perturbative 
Real-Space Renormalization Group (RSRG) |1^J, Con- 
tractor Renormalization Group (CORE) ^6|, or Dynam- 
ical Cluster Approximation (DCA) for embedded 4-site 
clusters coupled to an uncorrelated bath 17j. Similar at- 
tempts at studying small clusters (such as the ground 
state studies of weakly coupled Hubbard dimers and 
squares by Tsai and Kivelson ^3) have begun recently 
and the lessons learned here will be invaluable to such 
studies and useful to the condensed matter community 
in general. Above all, these small clusters can be synthe- 
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FIG. 1: Various possible configuration mixing of electrons 
(below liaff fifling) tliat can be found in an ensembfe of 4-site 
cfusters. Note tliat the mixing of configurations is brouglit 
about by tlie temperature. 



sized and utilized for understanding essential many-body 
physics at the mesoscopic level and hence are useful in 
their own right. 



II. MODEL AND METHODOLOGY 

The single orbital 4-site minimal Hubbard Hamilto- 
nian, 

H = ^(c+ c,+i, + H.c.) + uY,n^^n,i - /i^ 5^(1) 

ia i i 

with hopping t and on-site interaction U is the focus of 
this work. Periodic boundary conditions are used for 
the clusters. Unless otherwise stated, all the energies re- 
ported here are measured in units of t (i.e. t has been 
set to 1 in most of the equations that follow) . Our cal- 
culations are based on the exact analytical expressions 
for the eigenvalue En of the n*'' many-body eigenstate 
of the 4-site Hubbard clusters T^. As we show here, 
this model, used in conjunction with the grand canonical 
and canonical ensembles, yields valuable insight into the 
physics of strongly correlated electrons. 



A. Thermodynamics and response functions 

The complete phase diagram of interacting electrons 
can be obtained with high accuracy due to the ana- 
lytically (exact) given thermodynamic expressions. In 



3 



Fig-H possible electron configurations (below half- filling) 
in the grand canonical ensemble for the 4-site clusters 
are shown. The grand partition function T^u (where the 
number of particles iV and the projection of spin are 
allowed to fluctuate) and its derivatives are calculated 
(exactly) without taking the thermodynamic limit. The 
exact grand canonical potential for many-body inter- 
acting electrons is 



n<NH 



En-f^Nn-hs 



(2) 



where Nn and are the number of particles and the 
projection of spin in the n*'' state respectively. The 
Hilbert space dimension in ((SJ is Nh = 4^. The deriva- 
tives we study may be labeled as first order (such as 
the average spin projection/magnetization in response to 
an applied magnetic field) or second order (such as fluc- 
tuations/susceptibilities). These responses, evaluated as 
functions of chemical potential fi, applied field h, on- 
site Coulomb interaction U and temperature T, carry a 
wealth of information that can be used to identify vari- 
ous phases and phase boundaries. Some of these results 
for the 2- and 4-site clusters were reported earlier ll^ . 

The (first order) responses due to doping and external 
magnetic field are as follows: 



(N) 



9/i 



dflu 
'~dh 



(3) 



(4) 



Analytical expressions derived for the averages {N) and 
(s^) are analyzed numerically in a wide range of U , h, 
fi and T parameters. The charge and spin degrees re- 
spond to an applied magnetic field (h) as well as electron 
or hole doping levels (i.e. chemical potential fi) and dis- 
play clearly identifiable, prominent peaks, paving the way 
for rigorous definitions of Mott-Hubbard (MH), antifer- 
romagnet ic ( AF), spin pseudogaps and related crossover 
behavior . The exact expressions for charge suscepti- 
bility, Xc = and spin susceptibility, Xs = ^^T" '^^^ 
be found as a function of U, h, fi and T from, 



(X^) - {Xf = T 



dx 



(5) 



where X corresponds to TV or and x to or h. Using 
maxima and minima in spin and charge susceptibilities, 
phase diagrams m a T vs fi plane for any U and h can 
be constructed. This approach also allows us to obtain 
QCP and rigorous criteria for various transitions, such 
as the MH crossover at half-filling and MH like bifur- 
cations, using the evolution of peaks in charge or spin 
susceptibility ^| (see below). 



B. Charge (pseudo) gap 

We define canonical energies /i± , 

/i+ = E{M + 1, Af '; U :T)- E{M, M'; U, T) (6) 
^_ = E{M, M'; U -.T)- E{M - 1, M'; U : T) (7) 

where E{M, AI';U : T) is the canonical (ensemble) en- 
ergy with a given number of electrons N = M + M' 
determined by the number of up (M) and down (M') 
spins. At zero temperature the expressions and {Tj) 
are identical to those introduced in "2^. At finite tem- 
perature, the calculated charge susceptibility is a differ- 
entiable function of N and fi. The peaks (i.e. maxima) in 
Xc(m)i which may exist in a Umited range of temperature, 
are identified easily from the conditions, Xci^^±) — with 
Xc(/^±) < 0. We define Tc(/i) to be the temperature at 
which a (possible) peak is found in Xcip) at a given /i. 

The positive charge gap for electron-hole (excito) ex- 
citations, A''^''(r) > 0, is defined by 



A'=-''(r) 







if fi+ > fi- 
otherwise. 



(8) 



as the separation between and /i_ . The electron- 
hole instability A''-''(r) > can exist in a limited range 
of temperatures and U > Uc{T), with A'=~''(T) = at 
U < Uc{T). (In general, the critical parameter Uc{T), 
identified and discussed in Sections llll Al and 1111 ii\ de- 
pends also on h.) The energy gap, A'^^''(T) > 0, serves 
as a natural order parameter in a multidimensional pa- 
rameter space h,U,T and at T ^ will be called a pseu- 
dogap, since x has a small, but nonzero weight inside the 
gap. At T = 0, this gap A^j"'' = A'=-''(0) will be labeled 
a true gap since Xc is exactly zero inside. 

The difference — ^_ is somewhat similar to the dif- 
ference I— A for a cluster, where / is the ionization poten- 
tial and A the electron affinity. For a single isolated atom 
at half-filling and T = 0, I — A is equal to U and hence 
the above difference represents a screened U, reminiscent 
of Herring's definition of U in transition elements [2H. 



C. Mott-Hubbard crossover 

The thermodynamic quantities with fixed N {canonical 

approach) are certainly smooth, analytical functions of T 

and U. Thus, one may naively think that at half-filling 

and large [/, there is no real cooperative phenomena at 

T ^ U for transition from localized to delocalized elec- 
,2 

trons or at T jj for transition from antiferromagnetic 
to paramagnetic state. At finite temperature, the ther- 
modynamic quantities rtjj and (N) are both analytic and 
smooth functions of T and /i. Although the charge sus- 
ceptibility Xc is also a differentiable function of /i and T at 
all T > 0, Xc i>s n exhibits a weak, fourth order singular- 
ity at some critical temperature Tmh (saddle point) 1^. 
Thus the MH crossover at half-filling (/i — U/2) can be 
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defined simply as a critical temperature T]\jh, at which 
two peaks merge into one with /i — fi^ = /i_ = U/2 

and Xcif^) = ^^'^ Xc (a*) = 0; i-^- the tempera- 
ture corresponding to a point of inflexion in Xcil^) E3- 
This procedure gives a rigorous definition for the MH 
crossover temperature Tmn (from a localized into an itin- 
erant state), at which the electron- hole pseudogap melts, 
i.e. A'^-'^ITmh) = 0. The MH crossover, due to its 
many-body nature, is also a cooperative effect which may 
occur even for a single atom. 



D. Spin (pseudo) gap 

Analogously, we define a (positive) spin pairing gap 
between various spin configurations (projections of spin 
s) for a given number of electrons {N — M + M') in the 
absence of a magnetic field {h = 0) as, 

A%T) = E{M + 1, M' - 1; [/ : T) - E{M, M'; U : T).{9) 

This corresponds to the energy necessary to make an ex- 
citation by overturning a single spin. Possible peaks in 
the zero magnetic field spin susceptibility Xs{p)i when 
monitored as a function of /i, can also be used to define 
an associated temperature, Ts{p), as the temperature at 
which such a peak exists and the spin pseudogap as the 
separation (distance) between two such peaks. 



E. AF (pseudo) gap and onset of magnetization 

Similar to the (charge) plateaus seen in {N) vs ^, 
we can trace the variation of magnetization (s^) vs an 
applied magnetic field h and identify the spin plateau 
features, which can be associated with staggered mag- 
netization or short range antiferromagnetism. We cal- 
culate the critical magnetic field hc± for the onset of 
magnetization (s^ iO), which depends on N and /i, 
by flipping a down spin, hc+ = E{M + 1,M' — 1;U : 
T) - E{M, M'; C/ : T) or an up spinji^- = E{M, M'; U : 
T) - E{M - 1,M' + 1;U : T) The spin singlet 

binding energy A^^{T) > can be defined as, 



F. Pairing instability 

To determine whether the cluster can support electron 
pairing at finite temperature despite the purely repul- 
sive electronic interactions, the electron-electron (or hole- 
hole) pair binding energy, 

A^(r) = 

[E{M - 1, M'; U -.T)- E{M + 1, M'; U : T)] - 
2[E{M, M'; U -.T)- E{M + 1, M'; U : T)], (11) 

is calculated by adding or subtracting one electron near 
N ^ M + M'. The average energy E{M, M'; U : T) 
is given for configurations with a fixed number of elec- 
trons N and spin magnetization = using our grand 
canonical ensemble approach. At zero temperature, the 
binding energy (|ll|l is identical to the one introduced in 
Rcf. '23]. 

Using the definitions for /i± from Eqs. and |7| 
electron-electron (or hole-hole) pair energy can also be 
written as, 



A^(T) = 







otherwise. 



(12) 



In the ground state, the electron-pair binding energy 
gap at {N) sa 3 is fully developed at J7 < Uc{T) when 
A^(0) > 0, i.e. > /i-i., which leads to the phase 
separation into (N) = 2 and (N) = 4 clusters (see Sec- 
tion ^^Sl and an effective attraction between the elec- 
trons in (N) — 2 cluster configuration On the 
other hand, when fi^ > A*-; condition A^~^{T) > 
with U > Uc(T) provides an electron-hole pairing mech- 
anism as a precursor to antiferromagnetism p jj/ . We can 
define Tp as the temperature at which the pseudogap, 
A^(Tp) = 0, vanishes. The existence of particle-particle 
(A^(r) > 0) or particle-hole (A^-^(r) > 0) instability 
and the corresponding solution for a positive pseudogap 
(A(r) > 0) can be formulated at an arbitrary ?7 > by 
combining both equations ^ and (|12|l in one as. 



A(T) 



A^(T) ifC/<C/c(T). 



(13) 



.AF 



(T) 



hc+ — he 



(10) 



At zero temperature, A(0) = at C/ = or f7 = Uc{0). 



and serves as a natural antiferromagnetic order parame- 
ter in a multidimensional parameter space U, T, fi. This 
will be called an AF pseudogap at nonzero temperature. 
We define Tj^p as the temperature at which the pseu- 
dogap, A^-^{T) = 0, vanishes and above which a para- 
magnetic state is found. An exact analytical expression 
for the AF spin gap in the ground state (A^^(O)) at 
half- filling was obtained in Ref. In what follows, 

all of the temperatures defined above, Tc(/Lt), Ts(^) and 
Ta_f(m), will be used when constructing phase diagrams. 



III. RESULTS 
A. (A'^) and (s^) vs /i and pseudogaps 

In Fig. 13 we explicitly illustrate the variation of (N) < 
A vs ^ for various U values in order to track the variation 
of charge gaps with U . The opening of such gaps is a 
local correlation effect and clearly does not follow from 
long range order, as exemplified here. The true gaps at 
{N) = 1 and {N) ~ 4 develop for infinitesimal C/ > 
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average electron number and spin projection 
tT = 0.02 




-3.0 -2.0 -1.0 0.0 1.0 2.0 

chemical potential (i^) 



FIG. 2: Variation of average electron concentration (A'^) (top) 
and average spin (s^) (bottom) vs ji for various U values at 
temperature T = 0.02. The vanishing of the charge gap near 
(A'^) = 3 for U = A has implications related to pairing as 
discussed in the text. Note also that the spin plot has been 
obtained with an applied magnetic field of ft = 0.1 and shows 
the stabilization of a magnetic state near (A'') = 3 for (7 = 4. 
At zero field, (s^) = everywhere due to degeneracy between 
spin up and down states. 



and increase monotonically. In contrast, the charge gap 
at {N) = 3 opens at finite U > Uc{T) (see Fig. 121. Thus 
at low temperature, {N) (expressed as a function of fi in 
Fig- El evolves smoothly for U < Uc{T), showing finite 
leaps across the MH plateaus only at (TV) = 1, (TV) = 2. 
Such a density profile of (TV) vs ^ near {N) = 3 closely 
resembles the one calculated in Fig. |2 for the attractive 
4-site Hubbard cluster with [7 = -4 at T = 0.02 and is 
indicative of a possible particle pairing instability. At 
larger U > Uc{T), an electron- hole gap is opened at 
{N) — 3. Therefore the cluster at J7 > Uc{0) behaves 
as a MH like insulator at all (allowed) integer numbers 
(1 < TV < 8) with electron charge localized (non- Fermi 
liquid). In contrast, at C/ < Uc{Q) the chemical potential 
gets pinned upon doping in the midgap states at (N) ~ 
3. 

While Fig. |2 shows the magnetization at a relatively 
high magnetic field, its behavior at very low tempera- 
ture, r — > and infinitesimal magnetic field, ft, — > 
is also noteworthy. In this U increases, {N) 

and (s^) vs fj, near (TV) = 3 reveal islands of stability, 
due to phase separation (see Section UlI C|) . for various 
charge {N — 2 and iV = 4) and spin (s and s^) con- 
figurations as follows. Phase A {U < Uc{0)): particle- 
particle A^(0) > and spin-spin A^^'^^^O) > pair- 



ing gaps with the minimal projection of spin (s^) — 0. 
Phase B (£7^(0) < C/ < 4(2 + \/7) ~ 18.583): particle- 
hole A'=-''(0) > and spin-spin A^'=i/2(o) > pair- 
ings with the spin s = 1/2 and unsaturated ferromag- 
netism, (s^) = 1/2 {triplet pairing) Phase C (large 
U > 4(2-^^7)): particle-hole A'=-^(0) > pairing with- 
out spin gap (A*''=^/^(0) = 0) and maximum projection 
(s^) = 3/2 {saturated ferromagnetism) [25l |. 

In Phase A for U — 4, charge and spin are coupled 
{i.e. no charge-spin separation), while in Phase B at 
U — 6, the charge and spin are partially decoupled {par- 
tial charge-spin separation) . In Phase C for U —f oo, the 
charge and spin are fully decoupled, when the charge gap 
saturates to its maximum value, — > 2(2 — V2), while the 
spin gap from (s^) = 1/2 to (s^) = 3/2, defined earlier in 
0, vanishes {full charge-spin separation) . Phase A, due 
to strong particle-particle coupling with double electron 
charge {Q = 2e) and zero spin (s^ = 0) with a majority 
of {N) = 2 clusters, becomes a good candidate for the full 
'bosonization' of charge and spin degrees of freedom and 
possible 'superconductivity'. In contrast, at even num- 
bers of electrons, such as (N) ~ 2 in Fig. [3 there are 
electron-hole A''~''(0) > and spin-spin A'*"=°(0) > 
pairings at all U values, and therefore the charge Q — 2e 
and spin = are both coupled and there is full charge- 
spin reconciliation, when the singlet-triplet spin excita- 
tion gap at quarter filling approaches the charge gap, 
A^'=o = A"-'' = 2{2V2- 1), as f7 ^ oo. Exactly at half 
filling, (N) = 4, there is partial charge-spin separation at 
all finite U > Uc{0). However, as t7 ^ cx) the charge MH 
gap A^-''(0) oo and the AF gap A"'=0(0) (van- 
ishes); there is full charge-spin separation with saturated 
spin (s^) = 2 in this limiting case. Also, we find that for 
all U, the clusters with a single electron at (N) = 1 is a 
MH like insulator: charge gap A'^~''(0) oo, while the 
spin gap A*''=-^/^(0) with saturated spin (s^) — 1/2. 
For any given N in the charge sector, one can easily iden- 
tify an insulator or a metallic liquid if A'^~''(0) > or 
^e.-h(Q^ = respectively. Accordingly, it is also useful to 
distinguish a spin insulator, A''(0) > 0, or a spin liquid 
(A*(0) = 0) state in the spin sector. 

In Fig. we show the evolution of charge susceptibil- 
ity Xc as a function of /i, which exhibits clearly identifi- 
able sharp peaks. At low temperature, peak structures in 
Xc(m) a-nd zero magnetic field spin susceptibility, Xs{f'), 
are observed to develop in these clusters; between two 
consecutive peaks, there exists a pseudogap in charge or 
spin degrees. The opening of such distinct and separated 
pseudogap regions for spin and charge degrees of free- 
dom (at low temperature) is a signature of corresponding 
charge and spin separation away from half-filling. 



B. Pairing gap at (A'^) = 3 

At relatively large U > Uc{T), the energy gap A§(C/ : 
T) = E{A; U ■.T) + E{2; U : T) - 2E{3; U : T) becomes 
positive for {N) — 3 (see Fig. 0J. Its zero temperature 
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charge susceptibility peak positions at T = 0.02 energy gap vs U from Canonical ensemble 




FIG. 3: The charge susceptibility vs chemical potential jj. at 
T = 0.02 for two different U values ([/ = 4 and (7 = 6) below 
half-filling. Note that there are clearly identifiable peaks, and 
such peak positions, when monitored as a function of tem- 
perature, have been used to construct phase diagrams (see 
text). 



value was first derived analytically [13| as, 

AUU:T = 0) = --|V(16i2 + {/2)cos | + 
I _ 2^(48^2 + [/2)cos| + 
\/ 32*2 + [/2 + 4^64t4 + 3t^U^, (14) 

where a = arccos{(^^ - ^)/{^ + ^)^} and 7 = 
arccos{(4f7i2)/(l^ ^)i}. Due to (ground state) 
level crossings j23| , the exact expression (|14|l is valid only 
in a limited range of f7 < 4(2 + V7). The critical value, 
Uc{0) = 4.58399938, at which A^{Uc : T = 0) = 0, re- 
ported in Ref. serves as an estimation of the accu- 
racy of the gap value, which is slightly different from a 
value obtained in Ref. p^ . 

When A§(C/ : T) = EiA;U : T) + E{2-U : T) - 
2i?(3; U : T) becomes negative for U < Uc{T) as shown 
in Fig. 01 the (N) = 3 states become less (energetically) 
favorable when compared with (N) — 2 and (N) — 4 
states. This is a manifestation of electron binding where, 
despite bare electron repulsion, electron pairs experi- 
ence an attraction [Tg,, ,23, 24]. We have also observed 
a similar vanishing of charge gaps for negative U val- 
ues (see Fig. where there is an inherent electron- 
electron attraction, supporting the above statement. For 
U > Uc{T), the gap in the electron- hole channel is pos- 
itive (i.e. A§(T) > 0) which favors excitonic, electron- 
hole pairing, similar to MH gap at half-filling. 



FIG. 4: The energy gap given by Eq. 11411 at zero tem- 
perature and its evolution at several nonzero temperatures 
as a function of U obtained from the canonical ensemble. 
This is directly related to the electron-electron pairing gap 
{U < Uc{T)) and electron-hole charge gap (U > Uc{T)) as 
defined in the text. Note that charge pairing is unlikely to 
occur above temperature T^O.OS. As previously, all energies 
are measured in units of t, the hopping parameter. 



4-site phase diagram: U=4 




FIG. 5: Temperature T vs chemical potential /i phase dia- 
gram for the four-site cluster aX U = 4 and h — 0, below half 
filling {fi < jj,o). Regions I and II are paramagnetic phases 
and quit e similar to the ones found in the two-site cluster 
(Ref. [l3|)j showing strong charge-spin separation. Phase III 
is a MH antiferromagnetic phase (with zero spin). However, 
note the (new) line phase (labeled P) which consists of charge 
and spin fluctuations. This is a new feature seen in the 4-site 
cluster which suggests the existence of electron-electron pair- 
ing and phase separation into hole-poor and hole-rich regions 
at low temperature. For U = A, pairing occurs below the 
temperature, Tp=0.076. 
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4-site phase diagram: U=6 




chemical potential 



FIG. 6: Temperature T vs chemical potential /i phase diagram 
for the four site cluster at f/ = 6 and h = Q, below half-filling 
(m ^ Po)- Regions I, II and III are quite similar to the ones 
found for U=4 in Fig.|Kl again showing strong charge-spin sep- 
aration. However, a charge gap opens as a new bifurcation (I 
and II phases) which consists of charge and spin pseudogaps, 
(replacing the equilibrium line phase P in the previous figure, 
i.e. no electron-electron pairing here: see text for details). 
Temperature labels similar to those shown in Fig. |3 may be 
used here but have been suppressed for clarity. 



C. Phase separation 

It appears that the eanonieal approach yields an ade- 
quate estimation of a possible pair binding instability in 
the ensemble of small clusters at relatively low tempera- 
ture. The competition between attraction and repulsion 
under hole-doping can lead to a 'microscopic phase sepa- 
ration', which comprises of an inhomogeneous structure 
of competing and coexisting hole-rich {N) = 2 (hole-hole 
or electron-electron pairing), hole-poor {N) — 4 (AF at 
half filling) and magnetic {N) — 3 clusters. From Eqs. 
© and (I?!), it is apparent that the condition /i_ = /^_|_ 
at [/ = C/c(0) and {N) — 3 defines a lower bound for 
the existence of a phase separation boundary that distin- 
guishes a charge-spin separated region from a charge-spin 
coupled one. 

If we neglect every second hopping term in the two di- 
mensional square lattice, the system can be thought of 
as an ensemble of decoupled 4-site non-interacting clus- 
ters 0,01, as shown in Fig.^ New and important fea- 
tures appear if the number of electrons {N) — 3 and total 
magnetization (s^) = are kept fixed for the whole sys- 
tem of decoupled clusters, placed in the (particle) bath, 
by allowing the particle number on each separate cluster 
to fluctuate. One is tempted to think that due to symme- 
try, there is only a single hole on each cluster within the 
{N) = 3 in ensemble. However, this statement reflects a 
simple average only. Due to thermal and quantum fluc- 



tuations in the density of holes between the clusters (for 
U < Uc{T)), it is energetically more favorable to form 
pairs of holes. In this case, snapshots of the system at 
relatively low temperatures and at a critical value {fj,p in 
Fig. |3J) of the chemical potential would reveal equal prob- 
abilities of finding (only) clusters that are either hole- rich 
(N) = 2 or hole-poor (N) = 4. 

However, at higher temperatures when pairing coexists 
with magnetic spin fluctuations, there exists a small win- 
dow of parameters which brings some stability to (N) — 3 
clusters above Tp. Thus the crossover from full separa- 
tion (segregation) to the coexisting magnetic {{N) = 3) 
and hole-rich {{N) = 2) phases develops smoothly and 
depends on the degree of disorder, i.e. temperature. In 
section IIII Al it was shown how the changes in the pa- 
rameter U affect the spin (magnetization) by changing 
the cluster configuration. 

Phase separation of magnetic {{N) — 3) and paired 
{{N) — 2) states can also be triggered by increasing the 
magnetic field. The phase separation (i.e. segregation) 
into separate magnetic {N) — 3 and hole-rich (N) — 2 re- 
gions seen here at ^_ > closely resembles the phase 
separation effect detected recently in super-oxygenated 
La2-2;Sr2;Cu04+j^, with various Sr contents Thus 
our results are consistent with the experimental observa- 
tion that a small (applied) magnetic field mimics doping 
and promotes stability of the magnetic phase (N) = 3 
over the superconducting state (N) ~ 2 near optimal 
doping. 

In addition, our calculated probabilities (from the 
grand canonical ensemble) at low temperature show that, 
in the interval ^+ < fj, < jj,p, (N) = 2 clusters become 
the majority; i.e. we observe both, pairing and phase sep- 
aration at low temperature. Note that phase separation 
here has a different origin and occurs at relatively weak 
coupling U < Uc far from the level crossing regime, at 
which the spin gap vanishes. Therefore, the mechanism 
of phase separation here is different from the one found 
in Refs. |2a . |27| at large U limit due to the spin density 
fluctuations (/i+ < /i_) and spin saturation. Indeed, the 
four-site cluster at large U > 4(2-|-\/7) reveals ferromag- 
netism in accordance with the Nagaoka theorem . 



D. Phase diagrams 

In Figs. El and El phase diagrams for the 4-site clus- 
ter ([/ = 4 and U = 6) are shown (see also Ref ^2^). 
These have been constructed almost exclusively using the 
temperatures, Tc(/i), Ts(/i) and Taf{i^), defined previ- 
ously. We have identified the following phases in these 
diagrams: (I) and (II) are MH like paramagnetic phases 
with a charge pseudogap separated by a phase boundary 
where the spin susceptibility reaches a maximum, with 
A^-^(T) > 0, A^^ = 0; at finite temperature, phase I 
has a higher {N) compared to phase II; Phase (III) is a 
MH like antiferromagnetic insulator with bound charge 
and spin, when /^"-^{T) > 0, A''(r) > 0, A^^{T) > 0; 
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(P) is a line phase ioi U — A with a vanished charge gap at 
(N) — 3, now corresponding to the opening of a pairing 
gap (A^(T) > 0) in the electron-electron channel with 
A§([/ : T) < 0. We have also verified the well known 
fact that the low temperature behavior in the vicinity of 
half-filling, with charge and spin pseudogap phases coex- 
isting, represents an AF insulator JJ|J. However, away 
from half-fiUing, we find very intriguing behavior in ther- 
modynamical charge and spin degrees of freedom. 

In both phase diagrams, we find similar paramagnetic 
MH (I), (II) charge-spin separated phases in addition to 
the AF (III) phase where spin and charge are bounded. 
In Fig. spin-charge separation in phases (I) and (II) 
originates for relatively large U (= 6) in the underdoped 
regime. In contrast, Fig.Oshows the existence (at U ~ 4) 
of a line phase (with pairing) similar to U < case with 
electron pairing (A^(r) > 0), when the chemical poten- 
tial is pinned up on doping within the highly degenerate 
midgap states near (underdoped) 1/8 filling. 



E. Quantum critical points 

Among other interesting results, rich in variety, sharp 
transitions at QCP near (N) — 3 are found in the ground 
state at U > Uc{0) between phases with true charge and 
spin gaps; for infinitesimal T — > 0, these gaps are trans- 
formed into 'pseudogaps' with some nonzero weight be- 
tween peaks (or maxima) in susceptibilities monitored as 
a function of doping (i.e. /i) as well as h. In the limiting 
case Tcip-c) 0, the QCP doping, ^c, defines a sharp 
MH like (AF) like transition with diverging Xc E3- 
the QCP doping, /i^, with Ts{^s) 0, the zero spin 
susceptibility, Xs, also exhibits a maximum. 

In Fig. El the critical temperature Ts{fi) falls abruptly 
to zero at the QCP doping, /z^ (true only for U > Uc{Oj), 
implying ^ that the (spin) pseudogap can exist indepen- 
dently of possible particle pairing in Fig. |H| In contrast, 
for U < Uc{0) and low temperature, there is no charge- 
spin separation near (N) = 3. Therefore in Fig. [31 at 
U — A [U < Uc{Q)) we do not observe any QCP associ- 
ated with fis or fic close to (N) = 3. Instead, Fig.jslshows 
the existence of a line phase (with pairing) similar to the 
attractive U < case with a spin pseudogap, which ex- 
ists only at finite temperature Ts{fi) > 0, and electron 
pairing (Tp > T,), when the chemical potential is pinned 
up on doping within the highly degenerate midgap states 
near (underdoped) 1/8 filling. 

We have also seen that a reasonably strong magnetic 
field can bring about phase separation and has a dra- 
matic effect (mainly) on the QCP at fig, at which the 
spin pseudogap disappears. It is evident from our exact 
results that presence of QCP at zero temperature and 
critical crossovers at Tc{n), T's(/x) and Taf{p) temper- 
atures, gives strong support for the cooperative charac- 
ter of existing phase transitions and crossovers similar 
to those seen in large thermodynamic systems at finite 
temperatures |29ll30( . 



F. Charge-spin separation 



Charge-spin separation effect is fundamental for un- 
derstanding of the generic features common for small 
and large thermodynamic systems. We formulated ex- 
act criteria when the charge and spin excitations are 
decoupled at U > Uc{T). There is controversy regard- 
ing the nature of MH and AF transitions and relation 
between their consequent critical temperatures |l3j . In 
Figs. [SI and [^ the charge, decoupled from spin degree 
of freedom, condenses at temperatures below Tc while 
AF spin correlations near half-filling are seen to develop 
at lower temperatures Taf{ij-) < Tc{^) 0|. However, 
in the limited range close to /i > /ic, there is a reverse 
behavior, TApifJ-) > Tc{fi). Electrons were until recently 
thought to carry their charge and spin degrees of freedom 
equally; however accurate studies of thermodynamic re- 
sponse functions in nanoscale clusters show that in real 
materials, these two degrees of freedom are relatively in- 
dependent from one another and can condensate at dif- 
ferent doping levels /Zc, /is and transition temperatures 
T^ifi), Tap in) and Tsin) shown in Fig. [31 and Fig. [HI 

The charge-spin separation is an unusual behavior 
of electrons in some materials under certain conditions 
permitting the formation of two independent (bound) 
electron-electron or electron-hole pairs (quasiparticles) in 
charge sectors and spin singlet and triplet states in spin 
sectors. The spin quasiparticle only carries the spin de- 
gree of the electron but not the charge, while the charge 
quasiparicle has spin equal to zero but its electric charge 
equals either zero (electron-hole pair) or a charge of two 
electrons (electron-electron pair). We find that at large 
U > Uc{T), clusters with localized charge are favored 
over itinerant ones. 

As an important footnote, in the noninteracting U = 
case shown in Fig. \7\ the charge and spin peaks follow 
one another (in sharp contrast to the U = 4 and 6 cases). 
In regions I and II, positions of charge (as well as spin) 
maxima and minima coincide indicating that there is no 
charge-spin separation, even in the presence of a mag- 
netic field. In the entire range of /x, the charge and spin 
fluctuations directly follow one another without charge- 
spin separation. Our detailed analysis of the (responses 
such as) variation of electron concentration (N), zero 
magnetic field magnetization (s^) vs fj, and various fluc- 
tuations shows that there is no charge-spin separation 
and both, the spin and charge degrees, closely follow 
each other. Thus, at C/ = the spin and charge de- 
grees are strongly coupled to one another. On the other 
hand, the charge-spin separation effect at [/ 7^ led to 
rigorous deflnitions of Mott-Hubbard like, antiferromag- 
netic, spin pseudogaps, particle-particle pairings and re- 
lated crossovers. 
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4-site phase diagram: U=0 



ft 



charge max 
spin max 




chemical potential 



FIG. 7: The single particle or 'noninteracting' {U = 0) case, 
illustrating the positions of charge and spin susceptibility 
peaks in a T — ^ space for the 4-site cluster at fi < 1 (half- 
filling is at /io — 0). Note how the charge and spin peaks 
follow one another. Even in the presence of a nonzero mag- 
netic field there is no charge-spin separation. 



IV. CONCLUSION 

In summary, we have illustrated how to obtain phase 
diagrams and identify the presence of temperature driven 
crossovers, quantum phase transitions and charge-spin 
separation for any J7 > in the four-site Hubbard 
nanocluster as doping (or chemical potential) is varied. 
Specifically, our exact solution pointed out an important 
difference between the U = 4 ([/ < ?7c(0)) and U — & 
(U > Uc{0)) phase diagrams at finite temperature in 
the vicinity of hole doping « 1/8 which can be tied to 
possible electron-electron pairing due to overscreening of 
the repulsive interaction between electrons in the former. 



The resulting phase diagram with competing hole-rich 
{{N) = 2), hole-poor {{N) = 4) and magnetic ((A^) = 3) 
phases captures also the essential features of phase sep- 
aration in doped La2-a;Sr3;Cu04+j, [Tlj . Our analytical 
results near (N) « 3 strongly suggest that particle pair- 
ing can exist at [/ < Uc{T), while particle- hole binding 
is presumed to occur for U > Uc{T). It is also clear that 
short-range correlations alone are sufficient for pseudo- 
gaps to form in small and large clusters, which can be 
linked to the generic features of phase diagrams in tem- 
perature and doping effects seen in the HTSCs. The ex- 
act cluster solution shows how charge and spin (pseudo) 
gaps are formed at the microscopic level and their be- 
havior as a function of doping (i.e. chemical potential), 
magnetic field and temperature. The pseudogap forma- 
tion can also be associated with the condensation of spin 
and charge degrees of freedom below spin and charge 
crossover temperatures. In addition, our exact analyti- 
cal and numerical calculations provide important bench- 
marks for comparison with Monte Carlo, RSRG, DCA 
and other approximations. 

Finally, our results on the existence of QCP and 
crossover temperatures show the cooperative nature of 
phase transition phenomena in finite-size clusters similar 
to large thermodynamic systems. The small nanoclusters 
exhibit a pairing mechanism in a limited range of J7, 
and T and share very important intrinsic characteristics 
with the HTSCs apparently because in all these 'bad' 
metallic high Tc materials, local interactions play a key 
role. As charge and spin fluctuations are relevant to the 
charge and spin susceptibilities energy fluctuations 
are related to the specific heats and these new results 
for the 4-site and larger clusters, which provide further 
support to our picture developed here, will be reported 
elsewhere. 

One of us (ANK) thanks Steven Kivelson for interest 
and helpful discussion by communicating the results of 
their preprint (Ref. ^ISJ). This research was supported 
in part by the U.S. Department of Energy under Contract 
No. DE-AC02-98CH10886. 
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